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Mott transitions in the two-orbital Hubbard model with different bandwidths are investigated 
at finite temperatures. By means of the self-energy functional approach, we discuss the stability 
of the intermediate phase with one orbital localized and the other itinerant, which is caused by 
the orbital-selective Mott transition (OSMT). It is shown that the OSMT realizes two different 
coexistence regions at finite temperatures in accordance with the recent results of Liebsch. We 
further find that the particularly interesting behavior emerges around the special condition 
U — U' and J = 0, which includes a new type of the coexistence region with three distinct 
states. By systematically changing the Hund coupling, we establish the global phase diagram 
to elucidate the key role played by the Hund coupling on the Mott transitions. 
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Strongly correlated electron systems with some or- 
bitals have been investigated extensively. ' In partic- 
ular, substantial progress in theoretical understanding 
of Mott transitions of multi-orbital systems has been 
made by dynamical mean- field theory (DMFT) 3,4 calcu- 
lations. 5-29 Among them, the orbital-selective Mott tran- 
sition (OSMT) 30 in the multi-orbital system with differ- 
ent bandwidths has been one of the most active topics in 
this context. Typical materials are Ca2- x Sr x Ru04 31 and 
La n+ iNi n 3n +i, 32:33 where the OSMT is suggested to be 
realized by the chemical substitution and the change in 
the temperature. These experimental findings have stim- 
ulated theoretical investigations on this issue. 21-30, 34:35 
It has been clarified by extensive works that in the two- 
orbital system with different bandwidths at zero tem- 
perature, the OSMT occurs in general and merges to a 
single Mott transition in a certain restricted parameter 
region. 22 ' 24:25 Also, it has been realized that the Hund 
coupling is an important key parameter to understand 
the correct nature of the OSMT. 28,36 

In contrast to the detailed analysis at zero tem- 
perature, systematic discussions at finite tempera- 
tures 21, 23,27,28 are lacking. In particular, the stability 
of the intermediate metallic phase may be important 
to discuss the possibility of the OSMT in real mate- 
rials. It is thus desirable to investigate the OSMT in 
the multi-orbital model at finite temperatures. Quite re- 
cently, Liebsch have studied finite-temperature proper- 
ties of the OSMT by DMFT combined with the exact 
diagonalization method, 28 where the OSMT occurs even 
at finite temperatures. 

In this paper, we give the detailed investigation of the 
OSMT in the two-orbital Hubbard model at zero and fi- 
nite temperatures. We make use of the self-energy func- 
tional approach (SFA) proposed by Potthoff, 37 which is 
efficient to study finite-temperature properties, to deter- 
mine the detailed phase diagram. The present study pro- 
vides the results complementary to Liebsch's in a certain 
parameter regime, 28 and moreover establishes rich phase 
diagrams, which include a new type of the coexistence 
region with three distinct states. We further clarify how 



the magnitude of the Hund coupling controls the nature 
of phase diagrams at finite temperatures. 

We consider the two-orbital Hubbard model with dif- 
ferent bandwidths, which is given by the Hamiltonian, 
H = H a + J2 t W< with 



Ho = 
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where ct(ej ao -) is the creation (annihilation) operator 
of an electron at the ith site with spin a{~], [) and or- 
bital a(= 1,2), and rii aa = c\ a<y Ci aiy . Here, t a denotes 
the hopping integral for orbital a, /i the chemical po- 
tential, U(U') the intra-orbital (inter-orbital) Coulomb 
interaction, and J the Hund coupling including the spin- 
flip and pair-hopping terms. In the following, we impose 
the condition U = U' + 2 J, which results from rotational 
symmetry of degenerate orbitals. 

In order to discuss the Mott transitions at zero and fi- 
nite temperatures, we use here the self-energy functional 
approach SFA, 37 which is based on the Luttinger-Ward 
variational method. 38 This approach allows us to deal 
with finite-temperature properties of the multi-orbital 
system efficiently, 39 where standard DMFT with numer- 
ical techniques may encounter some difficulties in a prac- 
tical computation when the number of orbitals increases. 

In SFA, one makes use of the fact that the Luttinger- 
Ward functional does not depend on the detail of the 
Hamiltonian Ho as far as the interaction term H' is un- 
changed. 37 This enables us to introduce a proper ref- 
erence system with the same interaction term. One of 
the simplest models for the reference system is explicitly 
given by the following Hamiltonian, 7i re f = J2i T~L 
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where oS' (aS) creates (annihilates) an electron with u 
spin and a orbital, which is connected to the ith site in 
the original lattice. This approximation is regarded as 
a finite-temperature extension of the two-site DMFT 40 
in some respects. In this paper, we fix the parameters 
eoa = 0, e a = jj, and n = U/2 + U' — J/2 to discuss the 
zero and finite temperature properties at half filling. By 
choosing the parameters V a to minimize the grand po- 
tential, dQ/dV a = (a = 1, 2), we can find a proper ref- 
erence Hamiltonian, which approximately describes the 
original correlated system. 

We note that the hybridization V a between a given site 
on the original lattice and a site in the reference system 
may be regarded as the renormalized bandwidth of the 
a-orbital band. For example, for small V a , heavy quasi- 
particles are formed in the a band, and at V a = the sys- 
tem is driven to the Mott insulating state. By calculating 
hybridizations V\ and V2, we thus discuss the stability of 
the metallic state in our two-orbital Hubbard model. In 
the following, we focus on the system with W\ — 2.0 and 
W2 = 4.0 to discuss Mott transitions, where W a is the 
bandwidth of a-orbital for the semi-circular density of 
states (DOS), p a (x) = A/nW a ^l - (2x/W a ) 2 . 

We begin our discussions with the Mott transitions 
at zero temperature for a typical choice of the ratio 
of the parameters, V = 0.5U and J = 0.25U. Figure 
1 shows the contour plots of the grand potential and 
the corresponding entropy per site S/L at the station- 
ary point, where L is the number of sites. When the 
Coulomb interaction is small (U = 2.3), the stationary 
point, where the grand potential is minimized, is located 
around (Vi,V 2 ) ~ (0.14,0.55). This implies that the ef- 
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Fig. 1. The entropy per site S/L as a function of U with fixed 
ratios U' = 0.5U and J = 0.251/ at T = 0. Inset shows the 
contour plot of the grand potential for (a)J7 = 2.3, (b) 2.8 and 
(c) 3.6, where crosses denote the minima of the grand potential. 



fective bandwidth for each orbital state is finite, stabi- 
lizing the metallic (M) phase. Increasing the interaction, 
the stationary point of the system touches the V^-axis as 



shown in the inset (b): the effective bandwidth for the 
narrower band V\ vanishes while V2 still remains finite. 
Therefore, the OSMT occurs, where one band is insulat- 
ing and the other is still metallic. We refer to this phase as 
the OSM phase in the following. It is seen that the resid- 
ual entropy in the OSM phase is S/L = log 2, implying 
the formation of free localized s = 1/2 spins. Further 
increase of the interaction induces the Mott transition 
in the other band to the Mott insulating (MI) phase, 
where V\ = V2 = as shown in the inset (c). In this 
case, the localized spins for two bands form the triplet 
state due to the Hund coupling, so that the residual en- 
tropy S/L = log 3 in the region U > 3.4. Note that the 
residual entropy at zero temperature is an artifact of our 
simplified SFA approach. This pathological behavior can 
be somehow improved, e.g. by taking into account spin 
fluctuations more accurately. At zero temperature, we 
can find no other local minima in the grand potential 
when the interaction is varied. Therefore, we clearly ob- 
serve double second-order Mott transitions in the multi- 
orbital system, which indeed confirms the claim of Koga 
et al. 22 for the OSM transition obtained by means of 
DMFT combined with the exact diagonalization. 

We now move to the phase diagram at finite temper- 
atures. Since each transition at zero temperature is sim- 
ilar to that for the single-orbital Hubbard model, 3 the 
first-order transition may occur at finite temperatures 
around each critical point, as shown in ref. 28. In fact, we 
find that when the system belongs to the metallic phase 
close to the critical point, two stationary points repre- 
senting the M and OSM states appear simultaneously 
once the system is at finite temperatures (T = 0.002), 
as shown in Fig. 2. The double-well structure causes the 
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Fig. 2. (a) The contour plot of the grand potential Q(Vi,V2) at 
T = 0.002, U = 2.3 and J = 0.25U. Two minima coexist, which 
correspond to the M and OSM states, (b) The entropy vs the 
interaction U at T = 0.002, J = 0.25(7. 



first-order Mott transition with hysteresis at finite tem- 
peratures (Fig. 2 (b)). At T = 0.002, as U increases, 
the stationary point for the metallic state disappears 
around [/ c2 <~ 2.36, where the Mott transition occurs to 
the OSM phase. On the other hand, as U decreases, the 
OSM phase realized is stable down to U c \ ~ 2.24. The 
first-order transition point U c ~ 2.33 for T = 0.002 is de- 
termined by the crossing point of the two minima in the 
grand potential. The obtained phase diagram is shown 
in Fig. 3. It is found that the two distinct coexistence 
regions appear around U ~ 2.4 and U ~ 3.3. The phase 
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Fig. 3. U — T phase diagram at J = 0.2517. There are two co- 
existence phases, which are featured by the triangular-shaped 
regions: (left)the metallic phase and orbital-selective Mott phase 
coexist and (right) orbital-selective Mott phase and insulating 
phase coexist. 
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Fig. 4. (a) Entropy S/L and (b) Specific heat C/L as a function 
of U in the crossover region, J = 0.25(7. 

boundaries U&, U C 2 and U c merge at the critical temper- 
ature T c for each transition. Similar phase diagram was 
recently obtained by Liebsch by means of DMFT with 
the exact diagonalization method, who gave better es- 
timates of the critical temperatures. 28 Nevertheless, our 
SFA treatment elucidates further interesting properties 
such as the crossover behavior among the competing M, 
OSM and MI phases. To make this point more explicit, 
we calculate thermodynamic quantities such as the en- 
tropy and the specific heat, which are displayed in Fig. 4. 
We can see a double-step structure similar to that found 
at zero temperature in the curve of the entropy. Such 
anomalies are more clearly seen in the specific heat. It is 
quite impressive that the crossover behavior among three 
phases are clearly seen even at higher temperatures. We 
can thus say that the OSM phase is rather well defined 
even at higher temperatures above the critical tempera- 
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tures. 

We have so far discussed the OSMT for the system 
with rather large Hund couplings. According to the zero- 
temperature analysis, 22, 24, 25 it is known that around the 
special condition U = U' and J = at W\jWi = 0.5, 
a non-trivial single Mott transition occurs, whose origin 
is attributed to enhanced orbital fluctuations. 23 There- 
fore, it is interesting to observe what happens at finite 
temperatures, when the system approaches the condition 
U = U' and J = 0. By performing similar calculations, 
we end up with the phase diagram for the system with 
U' = 0.9417, J = 0.03C/ (Fig. 5). Since the Hund cou- 
pling is very small in this case, two critical points for 
the OSMT at zero temperature are close to each other. 
One of the most remarkable findings here is that three 
competing states (M, OSM and MI states) coexist in a 
certain region at finite temperatures, where three distinct 
local minima appear in the grand potential, as shown in 
the inset. In these parameters, the coexistence regions 




Fig. 5. U — T phase diagram at J = 0.03C7: we have used the 
same symbols as in Fig. 3. Note that the two coexistence re- 
gions arc not separated here in contrast to Fig. 3. Three dis- 
tinct states (M, OSM, and MI states) coexist in the region which 
is surrounded by filled and open circles. Inset shows the con- 
tour plot of the ground potential in the coexistence region for 
U = 4.45, T = 0.004, where there are three local minima. 

are quite sensitive to how large the Hund coupling is. 

To observe how the Hund coupling controls finite- 
temperature properties, we show the global phase dia- 
grams in Fig. 6 by systematically changing J. As the 
Hund coupling decreases, both the critical points 

(2) 

and U c2 increase monotonically, stabilizing the metallic 
state up to fairly large Coulomb interactions, as already 
mentioned for the zero-temperature case. 22 A remark- 
able feature in Fig. 6 is that the completely different 
behavior emerges in the critical temperatures and 
Tc 2 \ It is seen that T^p is almost unchanged, whereas 

(2) 

Tc strongly depends on the Hund coupling. This may 
be explained as follows. Around the first boundary be- 
tween the M and OSM phases, the nature of the tran- 
sition is essentially the same as the single-orbital case, 
as discussed above. Therefore, the Hund coupling lit- 
tle affects the critical temperature Tj 1 ^. On the other 
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Fig. 6. Overall picture for the coexistence regions with system- 
atic change of J: two types of the regions are shown for J = 
0.25(7, J = 0.1(7, J = 0.03(7 and J = from left to right. At 
J = 0, the transition from the metal to OSM disappears. 



hand, around the second phase transition between the 
OSM and MI phases, orbital fluctuations should play a 
vital role. The presence of the Hund coupling somewhat 
suppresses orbital fluctuations due to the formation of 
the triplet state. 15 Therefore, as J decreases, enhanced 

(2) 

orbital fluctuations raise the critical temperature Tc , 
as seen in Fig. 6. Note that for very small Hund cou- 
plings (J < 0.03[/) the OSM phase disappears, 25 and 
thus the finite as well as zero-temperature properties can 
be approximately described by the multi-orbital model 
(U + U') with the same bandwidths. 15 - 17 < 19 - 39 

In summary, we have investigated Mott transitions 
in the two-orbital Hubbard model with different band- 
widths at finite temperatures. By means of the self- 
energy functional approach, we have confirmed that two 
distinct coexistence regions appear in the phase dia- 
gram in accordance with the recent results obtained by 
DMFT. A notable new finding is that the OSM phase 
is rather well defined even at temperatures higher than 
the two critical temperatures. Further remarkable fea- 
tures have been elucidated in the phase diagram around 
special conditions with small Hund couplings J ~ 
(U ~ U'), where the coexistence region with three com- 
peting phases emerges. Orbital fluctuations are enhanced 
there, and therefore the system gets very sensitive to 
small perturbations, giving rise to the rich phase diagram 
at low temperatures. We have also clarified the role of the 
Hund coupling on global features in the phase diagram: 
the Hund coupling has little (sizable) effect on the tran- 
sition between the M and OSM (OSM and MI) phases, 
which is again attributed to the role played by enhanced 
orbital fluctuations. 

In this paper we have exploited the simplified version 
of SFA to elucidate fundamental properties at finite tem- 
peratures. Since instabilities to possible ordered phases 
have not been considered here, SFA is to be generalized 
to incorporate such ordered phases, which should be done 
in the future work. 
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